Disorder resistivity of solid neutron-star matter 
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Lower limits are found for the disorder electrical resistivity of solid neutron-star matter in the 
neutron-drip region which is amorphous and heterogeneous in nuclear charge. This temperature- 
independent resistivity, large compared with that produced by phonon scattering, has direct conse- 
quences for theories of neutron-star magnetic field generation and evolution. 
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Following the discovery of neutron stars in 1967, very 
many papers have been published containing theories of 
the nature and origin of their magnetic flux distributions. 
The original view was that the distributions are fossil 
fields amplified during gravitational collapse. A recent 
variant of this idea is field generation by rapid dynamo 
action during a brief interval of convective motion im- 
mediately after that process 0. But some authors have 
doubted that amplified fossil fields can be of the order 
of magnitude inferred from measured pulsar periods and 
their time-derivatives, and have considered thermoelec- 
tric field generation in the solid crust and liquid surface 
region of the star |^ |E Q during the early years after 
its formation. This Letter shows that disorder electrical 
resistivity is present in solid neutron-star matter. It is 
temperature-independent and so large that, with the im- 
purity resistivity, it severely constrains other than fossil- 
field theories. 

At times when thermoelectric processes have become 
negligible, and with no movement of matter, the evolu- 
tion of the magnetic flux density B in the solid crust is 
given by the Maxwell induction equation with Ohm's law 
and neglect of displacement currents. 
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where a is the zero-field conductivity and n the electron 
number density. Eq. (1) assumes that B is not so large 
in relation to n that classical conductivity fails (see, for 
example, Ref. 0). The existence of disorder resistivity 
leads to small values of a, renders negligible the term in 
(1) bilinear in B representing Hall drift, and allows un- 
ambiguous calculation of field evolution within the crust. 

Previously, values of a have been derived from calcu- 
lations of electron-phonon scatteriiig in a homogeneous 
bcc lattice , and from scattering 0, S by fractional 
concentrations Ci of impurity charges with mean Z 
and mean square deviation Q = J^i^i {Zi ~ Z) . The 



latter has been treated as a free parameter, usually as- 
sumed to satisfy Q <^ 1. For these cr, the Hall term 
in (1) is, in man y cases, larger than the ohmic dissipa- 
tive linear term [lOl Eq. (1) has some similarity 

with the vorticity equation for a viscous incompressible 
fluid (see, for example, Ref. 0) and a number of au- 
thors |ll[ll,l[ll4l | have investigated the existence of a cas- 
cade to higher-wavenumber components analogous with 
the Kolmogorov cascade in fluid turbulence. This would 
have a signiflcant effect on evolution because, although 
Hall drift is not in itself dissipative, the ohmic dissipation 
rate for any flux-distribution component is dependent on 
the square of its wavenumber. 

It must be emphasized that the references cited so far 
are not complete but provide an introduction to a very 
large literature on a subject of some complexity. Solu- 
tions of Eq. (1) are not simple if the Hall term is ap- 
preciable. Also, if a is derived from electron-phonon re- 
sistivity, it is so highly temperature-dependent that the 
cooling of the star and the evolution of its magnetic flux 
distribution must, in principle, be calculated in parallel. 
The present work shows that disorder becomes the dom- 
inant source of resistivity as the interior cools below 10^ 
K and so obviates these problems. 

There is now evidence that the assumption of a bcc lat- 
tice, on which calculations of a have been based, is not 
correct. Weak- interaction equilibrium is not reached dur- 
ing the formation and cooling of the neutron star so that, 
as a consequence of shell structure and of the small differ- 
ences in formation enthalpy between nuclei of different Z, 
the solid is both heterogeneous and amorphous [islll^. 
(This work applies to neutron stars that have not under- 
gone signiflcant mass-accretion. The high mass-transfer 
rates in certain binary systems can completely replace the 
natural solid crust with accreted matter of heterogeneous 
Z 17J.) It is well known that amorphous systems have 
high and substantially temperature- independent resistiv- 
ities (see, for example, Ref. ^3). In solid neutron-star 
matter, disorder is the most important source of resistiv- 
ity and does not seem to have been considered previously. 
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If the weak-interaction strength were (hypothetically) 
adjusted so as to allow the attainment of equilibrium, 
an amorphous and Z-heterogeneous solid would trans- 
form to a bcc structure with homogeneous nuclear charge 
and a population of lattice defects. Thus an amorphous 
structure with minimum disorder can be obtained by 
the reverse transformation, starting from a defect-free 
bcc lattice and sequentially changing each nuclear charge 
Z ^ Zi randomly so as to give fractional concentrations 
Ci. This procedure has been carried out for a system of 
24357 bcc lattice sites contained within a spherical vol- 
ume using distributions Ci of Zi found by the methods 
ofRef. |llE3- 

These are based on equilibrium bcc lat- 
tices with charge Zqldm whose parameters have been 
obtained at temperature T = for each of the matter 
densities in Table I using the compressible liquid-drop 
model of Lattimer et al [l9j . The densities selected are in 
the neutron-drip region and represent an interval equiv- 
alent to most of the depth of the crust. Calculations of 
formation enthalpies for nuclei with 20 < Zj < 50 in equi- 
librium with the neutron continuum |15[I1(tI | give values of 
the Ci for thermal equilibrium, at an initial temperature 
To — 5 X 10^ K, which is near the melting temperature 
Tm of the solid. There is movement away from weak- 
interaction equilibrium during the early stages of cooling 
and the fractional concentrations Ci evolve to frozen-in 
values giving the parameters Q of Table I. This produces 
a small shift 6Z in nuclear charge so that the mean charge 
of Table I is Z ~ Zcldm + SZ. Because there are con- 
siderable uncertainties in obtaining the shell-corrected 
nuclear formation enthalpies, a range of possibilities is 
considered Rows labelled sp in Table I are based 

on values of the Ci found from formation enthalpies given 
by the Strutins ki p rocedure for shell-correction (see, for 
example, Ref. l20|) with proton pairing. Rows labelled 
p are for Ci given by formation enthalpies with proton 
pairing but no shell correction. 

Each change of nuclear charge Z ^ Zi produces dis- 
placements of all other lattice sites in the system con- 
sidered. Our assumption here is that these are identical 
with those found in Ref. for an isolated point-defect 
in a homogeneous bcc lattice. (This, of course, is not 
strictly correct because as the sequence of transforma- 
tions Z —>■ Zi is made, the neighbours of a given nucleus 
become heterogeneous in charge and their positions de- 
viate from those of the bcc lattice.) These displacements 
were calculated for each of the neighbour sets a = 1 — 3, 
but for a > 3, they are assumed equal to those produced 
by a point defect in an isotropic incompressible elastic 
medium. The volume change represented by these dis- 
placements is determined, principally, by the properties 
of the Coulomb-electron stress tensor in the neutron-drip 
region whose isotropic components are between one and 
two orders of magnitude larger than the off-diagonal, and 
by the fact that Ata ~ 1, where k is the electron screening 
wavenumber [isl . Thus the displacements and vol- 



ume change are such that the electron density, averaged 
over a volume of the order of a Wigner-Seitz cell, adjusts 
to values almost exactly equal to those of the undisturbed 
lattice with Fermi wavenumber kpg. The result of calcu- 
lating the total displacement for each lattice site in the 
system considered is that the radial density distribution 
of nuclear positions relative to the origin site changes 
from a sum of (5-functions in the undisturbed bcc system 
to a sum of gaussian functions of approximately constant 
variance w^. There is also a change in the mean spatial 
density of sites, but it is statistically consistent with zero. 
Values of w have been calculated for each matter density 
and are given in Table I. 

Our calculation of disorder resistivity follows closely 
the classic paper of Ziman jij] on liquid metals. Com- 
parison with the exact electron scattering matrix element 
|22| as a function of the momentum transfer q indicates 
that for the size of momentum transfer allowed in the 
neutron-drip region, < q < 2kpe, the Born approxima- 
tion is adequate for scattering by the potential inside an 
electrically neutral Wigner-Seitz cell. The Fourier trans- 
form of the charge density is ZeF{q). It assumes the 
uniform proton density of a CLDM nucleus |l9j | and uni- 
form electron density inside a charge-neutral spherical 
cell whose radius is defined by Z and kpe. The scat- 
tering amplitude for the whole system is then a sum of 
two terms; one dependent on Z and the other on the dif- 
ferences (Zi — Z). The static structure function derived 
from the first term. 
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which vanishes in the absence of disorder, has been ob- 
tained directly, to terms of order g^w^, from the Fourier 
transform of the radial density distribution of the nuclear 
positions Ts Nuclei of any specific charge Zi are randomly 
distributed and the set of such nuclei has static structure 
function Sz — 

There is no interference between the two amplitude 
terms and so the resistivity TZ is 



n = n, + nd = ^ — (qa, + z^a^) 

ZCfle 



(3) 



where He is the electron chemical potential. The modified 
Coulomb integrals are. 



dq 

<7 V 4fc 



F^{q)SdAq)^' 



-2W 



(4) 



in which the final term, the Debye- Waller factor, can be 
set equal to unity at all q < 2kpf. for the relevant T <C 
Tp) , where Tp, is the Debye temperature of a neutron-drip 
bcc lattice. The first term in Eq. (3) is identical with the 
standard expression for impurity resistivity but the Table 
II values of TZz obtained from it are smaller than those 
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TABLE I: Disordered systems at matter density p have mean 
nuclear charge Z and electron Fermi wavenumber kpe giving 
an equivalent bcc lattice constant a. The parameter Q is the 
mean square nuclear charge deviation and is the variance 
of the gaussian radial distributions of nuclear positions rel- 
ative to the origin site. Rows labelled sp are for the cases 
in which nuclear concentrations d have been obtained by the 
procedures of Ref. with shell corrections and proton pair- 
ing. Rows labelled p are for Ci given by proton pairing alone. 
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of Ref. |3, |2( owing to the different assumption in the 
present work concerning F{q) . (The assumption of 8], 
that F{q) is the nuclear form factor divided by the smaU-g 
Hmit of the longitudinal static dielectric constant, gives 
the bracketed TZ^ values.) The second is the disorder 
resistivity and, even for the minimum disorder case, it is 
the larger term given the F{q) of this Letter, as shown 
in Table II. 

It is difficult to be specific about examples of greater 
disorder than the minimum defined here. One possibil- 
ity, which we arbitrarily define as maximum disorder, 
would resemble the inital stages of component separation 
in which small volumes of homogeneous lattice exist as 
irregular polyhedra, with random orientations. It is as- 
sumed that their linear dimensions are so small that the 
wave function of an electron is almost independent of the 
reciprocal lattice vectors of the polyhedron within which 
it is moving. In this case, and under the assumption of 
random orientation, the structure function (2) becomes 

a 

where Sa is the structure function for a single polyhe- 
dron of Na nuclei. Order of magnitude estimates of the 
resistivity can be found by substituting the mean Sa for 
Sd in Eq. (4). Values of Sa have been calculated for 
a number of trial irregular polyhedra containing 10'^ 
nuclei. Typical values give the resistivities 7^™°^ in the 
final column of Table II, which are not sensitive to poly- 
hedron size. Their order of magnitude shows that more 
disordered structures can produce very high resistivity. 
On this basis, we regard the values of TZ in the penulti- 
mate column of Table II as lower limits, though it has 
to be emphasized that they are subject to considerable 
uncertainties arising from the shell-correction procedure, 
as is shown by the differences between those labelled p 
and sp. 



TABLE II: Resistivities are listed for systems with the prop- 
erties given in Table I, 7?.™*" from the effect of minimum dis- 
order and TZz from charge heterogeneity; TZ — TZ™" + TZz- 
The bracketed values of TZz are for the form factor assumed 
in The final column gives an estimate of the resistivity 
arising from maximum disorder as defined here. 
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sumption about electron wave functions would fail and 
the solid would become ordered over those length scales 
which determine resistivity. However, neutron star crusts 
differ from white dwarf interiors in which component sep- 
aration is possible (see, for example Ref. j23) because 
cooling to T w O.STm occurs rapidly, in 10^~'^s. Vacancy 
mobility (or its analogue in an amorphous system) is the 
dominant diffusion mechanism enabling component sep- 
aration to occur {24] . Under the assumption that suf- 
ficient sinks exist in the solid structure to maintain a 
thermal equilibrium vacancy concentration during cool- 
ing, the jump frequency for a given nucleus is of the order 
of VD e'xp{— P{H Fv + i)) , where is the Debye frequency 
and l3~^ = fc^T. In a bcc lattice, the jump activation en- 
ergy given in Ref. |2^ can be expressed, approximately, 
as e « 15^VvKS^^/4, in which /i is the shear modulus, 
Vws is the mean volume per nucleus and 5^ = 0.08. 
Evaluation of this expression, as an estimate of e for an 
amorphous system, gives values much smaller than cal- 
culated vacancy formation enthalpies Hp^ 25] which are 
therefore the principal factor determining diffusion rates. 
For these formation enthalpies, jump frequencies become 
negligible at T < Tm relative to cooling rates. A further 
consideration is the formation of clusters of like nuclei 
which might form the initial stages of component sepa- 
ration. Approximate calculations [iSl ] of the formation 
enthalpy of some specific forms of cluster show that any 
decrease in enthalpy attained is smaller than the free en- 
ergy term derived from configurational entropy at Tm- 
On these grounds, we believe that component separation 
in neutron-star matter is insignificant and that the amor- 
phous nature of the solid persists during cooling. 

The consequences of disorder resistivity are unambigu- 
ous. FromEq. (1), the exponential decay time for a mode 
of wavelength A is r = ctA^/ttc^. Under the assumption 
that a is derived entirely from the high-density TZ of Ta- 
ble II, the decay time is r = Tq = 0.7 — 1.3 x 10^ yr for a 
mode with A/2 = 10^ cm, the approximate depth of the 
whole crust in a 1.4AfQ neutron star. This limits the ap- 
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plicability of models in which the field is confined to the 
crust j26| or has been submerged within it by post core- 
collapse accretion IJ^. Non- fossil field models, in which 
considerable amplification occurs after formation of the 
neutron star, rely on processes which are optimized at 
small electron density n. This is true of thermoelectric 
amplification |^ Q| and of amplification arising from 

the Hall term in Eq^ (1) which requires the inequality 

aB ;» nec 0, 0, Q|. Even if such processes occur 
for the small n at depths ~ 10^ cm, the models rely on 
movement of flux by Hall drift to regions of higher n in 
which ohmic dissipation is very slow. Previous assump- 
tions about a were consistent with regions of this kind 
occupying most of the crust volume. But the magnitudes 
of resistivity produced by disorder give ohmic dissipation 
times which are small everywhere in the spherical nuclear 
phase of the crust Q . 

For different limiting regions of B, solutions of Eq. (1) 
represent either Hall drift or ohmic dissipation. A re- 
cent paper |29l |. which also summarizes earlier work, has 
defined the boundary between these regions as a func- 
tion of Q. The conclusion was that, for Q ^ 1 and 
B - 1013 G, HaU drift never dominates the evolution 
of currents in the crust. Given these authors' assumed 
value a « 4.4 x lO^'^Q"^ s"^, this is consistent with the 
order of magnitude condition aB < nec being applied 
in the high-density regions of the crust. This Letter has 
shown that for amorphous structures, resistivity is not 
defined directly by Q but rather by the disorder which is 
an inevitable consequence of Z-heterogeneity. The spe- 
cific values of TZ calculated here allow one to reach, but 
with no free parameter, the same conclusion as |^ . Thus 
field evolution in neutron stars with B < 10^^ G, includ- 
ing almost the whole population of radio pulsars, is gov- 
erned by movement of flux from the liquid core, if it oc- 
curs, and by ohmic dissipation in the crust. The Maxwell 
stress is at least of the order of the critical shear stress 
of the crust at depths ^ cm. Therefore Eq. (1) may 
not be satisfled if, during the evolution of B, there are 
departures from hydromagnetic equilibrium sufficiently 
large to cause mechanical failure of the solid. But this 
complication does not affect the main conclusion of the 
present Letter, which is that in neutron stars with ages 
exceeding a small multiple of Tq, the magnetic flux den- 
sity observed is a fossil field, irrotational in the crust. 

The fields of 10^^"^'^ G thought to be present in the 
anomalous X-ray pulsars (AXP) may evolve by Hall drift. 
There are several possible mechanisms for the persistent 
emission of these objects j^Ol involving departures from 
hydromagnetic equilibrium in the star. But the resistivi- 
ties of Table H mean that, unless the fleld is almost irro- 
tational in the solid, ohmic dissipation must be a source 
of strong thermal emission with a lifetime ~ 10^~^ yr 
whether or not hydromagnetic equilibrium exists. 
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